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Abstract 

We present an alternative derivation of upper-bounds for the bulk modulus of both two-dimensional and three-dimensional 
cellular materials. For two-dimensional materials, we recover exactly the expression of the Hashin-Shtrikman (HS) upper-bound 
in the low-density limit, while for three-dimensional materials we even improve the HS bound. Furthermore, we establish 
necessary and sufficient conditions on the cellular structure for maximizing the bulk modulus, for a given solid volume fraction. 
These conditions are found to be exactly those under which the electrical (or thermal) conductivity of the material reaches 
its maximal value as well. These results provide a set of straightforward criteria allowing to address the design of optimized 
cellular materials, and shed light on recent studies of structures with both maximal bulk modulus and maximal conductivity. 
Finally, we discuss the compatibility of the criteria presented here with the geometrical constraints caused by minimization of 
surface energy in a real foam. 



Cellular solids appear widely in nature and are man- 
ufactured on a large scale by man. Examples include 
wood, cancellous bone, cork, foams for insulation and 
packaging, or sandwich panels in aircraft. Material den- 
sity, or solid volume fraction, <f>, is a predominant pa- 
rameter for the mechanical properties of cellular materi- 
als. Various theoretical studies on the mechanical proper- 
ties of both two-dimensional (2D) and three-dimensional 
(3D) structures have been attempted Q. Unfortunately, 
exact calculations can be achieved for cellular materials 
with simple geometry only , and numerical simulations 
or semi-empirical models 4] 5] 6] are required in or- 
der to study the mechanical properties of more complex 
structures. However, expression of bounds on the effec- 
tive moduli can be established. Perhaps the most famous 
bounds are those given by Z. Hashin and S. Shtrikman for 
isotropic heterogeneous media 0@. In particular, the 
Hashin-Shtrikman bounds for the effective bulk modulus 
in the low-density asymptotic limit (<f> -C 1) read: 
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for 3D structures 0- k^ 2D ^ and k^ 3D ^ are the actual 
bulk modulus respectively for 2D and 3D structures, and 
E, G, K are the Young modulus, shear modulus, and bulk 
modulus of the solid phase, respectively. These three 
elastic moduli are related by: E = g§ for 2D bodies 
and by: E = for 3D bodies. 

The search for optimal structures maximizing some 
specific modulus (for a given value of solid volume frac- 
tion 0), is of evident practical importance. In a re- 
cent study, Torquato et al. identified values of 



conductivity and elastic moduli of the two-dimensional 
square, hexagonal, kagome and triangular cellular struc- 
tures, and observed that the bulk modulus of these struc- 
tures is equal to the HS upper-bound value. The au- 
thors did not attempt to explain this result, although 
they noticed that such structures under uniform com- 
pression deform without bend (affine compression). Are 
these structures the only structures with maximal bulk 
modulus ? And if they are not, can we provide criteria 
on the structure of " optimized" cellular materials ? More 
intriguingly, Torquato et al. noticed that these structures 
present maximal electrical (or thermal) conductivity as 
well. Is this feature caused by the particular symmetry of 
the studied strutures, or is there an underlying relation 
between the conductivity and the bulk modulus of cel- 
lular materials ? We shall answer to all these questions 
in the present study. Indeed, Durand & Weaire [ll||T^ 
already established necessary and sufficient conditions on 
the structure of cellular networks having maximal aver- 
age conductivity. A quite similar approach is used in this 
paper to show that the very same conditions are also 
necessary and sufficient to maximize the bulk modulus 
of an open-cell material. There are some evident sim- 
ilarities between the constitutive laws (Ohm's law and 
Kirchhoff's laws) of electrical current in wires and those 
of the thin beam theory, but complexity is increased in 
the latter case, the scalar quantities I (the electrical cur- 
rent) and V (the electrical potential) being replaced by 
the vectorial quantities F (force acting on a beam) and 
u (the displacement field). 

We consider first the case of a 2D cellular material. We 
suppose its solid volume fraction (f> is sufficiently low, so 
the cell edges can be approximated as thin beams. Beams 
a priori can be naturally curved and have non-uniform 
cross-sections, as long as the cross-sectional area (I) of 
each beam (i,j) (where I refers to the curvilinear coordi- 
nate along the beam, and i and j denote the two nodes 
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linked by the beam) is small compared with its length 
squared. Let us isolate a circular portion of this material, 
of radius i?, and impose a uniform radial displacement 
—5Re r on its boundary (e r is the radially oriented unit 
vector; the body is under uniform tension when 5R < 0, 
and under uniform compression when SR > 0, see Fig. 
0). We define the 2D bulk modulus k^ 2D ^ of this structure 
as: 
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where SP is the average load applied on the boundary. 
We must point out that the definition above is different 
from the usual definition of bulk modulus: in the usual 
definition, a uniform radial load is applied on the surface 
of the material, while in the present definition a uniform 
displacement is imposed on its surface (allowing to ex- 
tend the notion of bulk-modulus to non-isotropic mate- 
rials). However, the two definitions are identical for 2D 
materials with square or hexagonal symmetry and for 3D 
materials with cubic or isotropic symmetry. 




FIG. 1: Circular portion of a 2D cellular network subjected to 
a uniform radially oriented displacement —8Re r of its bound- 
ary. The network is made of thin beams with a priori nat- 
ural curvatures and non-uniform cross-sections. We associate 
a bulk modulus for such a strain, defined as: k' 215 - 1 = -jf|j> 
where SP is the average load applied on the boundary. 



The expression of an upper bound can be easily es- 
tablished using the principle of minimum potential en- 
ergy: among all kinematically admissible displacement 
fields ( i. e. any displacement field twice continuously dif- 
ferentiable satisfying the displacement constraints on the 
boundary), the actual displacement (i.e. the one sat- 
isfying the equations of mechanical equilibrium) is the 
one that makes the potential energy an absolute mini- 
mum. Let u* (r) be the displacement field which satis- 
fies the equations of equilibrium throughout the body 
and the conditions on the boundary, u (r) any kine- 
matically admissible displacement field, and U ({u* (r)}) 



and U ({u(r)}) the respective potential energy associ- 
ated with these two displacement fields. Then, according 
to the principle of minimum potential energy: 



U({u* (r)})<f/({u(r)}) 



(4) 



Let us choose as kinematically admissible displacement 
field: u (r) = — ^r, and let us evaluate the potential 
energy associated with. We assume the cross-section of 
each beam is sufficiently small so u (r) is uniform on it 
(or equivalcntly, we suppose u (r) is a macroscopic field 
which has a uniform value on the beam cross-section). 
Thus, the stress tensor expressed in the local orthog- 
onal coordinate system has only one non-zero compo- 
nent: the axial-axial component. Consider an infinitesi- 
mal piece of a given beam (i, j), of length dl and cross- 
sectional area Sij{l). We denote ym and YM+dr the 
position of its two ends. Their relative displacement 
(u (rM+dr) — u (tm)) is colinear to the local tangent 
unit vector ty = meaning that the piece of beam 
deforms by axial compression only. The force Fij(l) act- 
ing on the surface Sij(l) is parallel to tjj and given by: 
Fy(Z) = Esij(l)^tij, where E is the Young modulus of 
the solid phase. The strain energy associated with such a 
deformation is -j s ij(0 (if) ■ Invoking additivity of the 
potential energy : 



t/({u(r)}) 



dl 



(5) 



(where the discrete sum is carried out on all the beams 
(i,j)), and introducing the volume fraction of solid: 



we obtain: 



0,j)o 



U({u(r)}) = ^Ecf > (SR) 2 . 



(6) 



(7) 



On the other hand, the actual potential energy 
U ({u* (r)}) is equal to half the work done by the ex- 
ternal forces (l3| : 

U ({u* (r)}) = ^SP2ttRSR = 2nn {2D) (SR) 2 . (8) 

Comparison of Eqs. [7| and [S] finally leads to an upper- 
bound for the bulk modulus: 



< 



E(j> 



(9) 



The same argumentation can be used for 3D open-cell 
structures: in that case, we find that the bulk modulus 



J3D) 



associated with a spherical portion of material of 



radius R, and defined as: 

1 1 4ttR 2 SR _ 3 SR 

k( 3D ) ~ #7ri?3 SP ~ RSP' 



(10) 
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is bounded as it follows: 

< f . (11) 

The upper-bound value we obtain is then lower than the 
HS upper-bound valued giving a sharpest estimation of 
the actual bulk modulus value. 



Criteria for maximal bulk modulus 

The upper bounds established above are optimal 
bounds, i.e. there exist cellular materials with maximal 
bulk modulus value. Criteria on the structure of such 
materials can be provided. Indeed, we show in the fol- 
lowing that the bulk modulus equals the upper-bound 
value if and only if the three following conditions are 
simultaneously satisfied: 

a) All the edges are straight. 

b) Each edge has a uniform cross-section area: Sij (I) = 

Sij . 

c) Every junction (i) between edges satisfies 
J2j s ij e ij = 0) where are outward-pointing unit 
vectors in the directions of adjoining edges. 

The demonstration is straightforward: according to 
the principle of minimal potential energy and the unique- 
ness of the actual displacement field, the inequality 0] 
becomes a strict equality if and only if the trial displace- 
ment field u (r) = — is the actual displacement field 
satisfying the equations of mechanical equilibrium. In- 
spection of force and moment balances along each beam 
and at each junction leads to the three necessary and suf- 
ficient conditions stated above. Let us make this precise. 
Consider a specific beam (i,j) (see Fig. EJl; at equilib- 
rium, the moments of forces acting on it must balance. 
Choosing as referencing point for the moments the node 
i, and denoting r^/ = tm — r$, where and tm are 
the respective position vectors of node i and of any point 
M belonging to the beam, we obtain: r,^ x Fjj(Z) = 
Esij(l)^YiM xtij = 0. Thus, the tangent unit vector t,j 
must be parallel to the position vector r;M for any point 
M belonging to the beam, leading to condition a). The 
forces acting on any piece Al of the straight beam 
must balance as well: E^Sij(l) — E^Sij(l + Al) what 
immediately leads to condition b). Finally, mechanical 
equilibrium at every junction i is satisfied if: Y]j F,y = 0, 
with Fij = —Esij^eij, leading to condition c). The mo- 
ment balance at every junction is automatically satisfied 
when conditions a), b), c) are fulfilled, since the force 
acting on each straight beam is then axially oriented. 
Furthermore, we check that the geometrical constraint 
on angles between adjoining edges is also satisfied; if e- • 




FIG. 2: Schematic of a paticular beam I denotes the 

curvilinear coordinate of a given point M along the beam. 
TiM and ty are respectively the position vector of point M 
taken from node i and the local tangent unit vector. 

denotes the unit vector parallel to the beam (i,j) after 
deformation and = Vj — r$, then: 

r i3 -+u(r 3 -)-u( ri ) _ (1 - f ) Tjj 
« ||r« + u(r J 0-u(r < )[| (l - f ) ||r«|| 

(12) 

what proves the material deforms by affine compression, 
and the angles between edges are preserved. 

Some comments 

We first summarize the limits of the theory: the solid 
volume fraction is supposed to be low enough for the thin 
beam theory to be valid. Moreover the relative imposed 
displacement ^ must be small enough so that Hooke's 
law can be used and no mechanical Euler instability oc- 
curs when the body is under compression. 

Conditions a), b), c) are the necessary and sufficient 
conditions to maximize the average conductivity of a net- 
work of thin wires as well 0] . Why structures satisfying 
these conditions have both maximal bulk modulus and 
maximal conductivity ? When the three conditions are 
fulfilled, the force acting to each beam is then parallel 
to it, and the corresponding deformation of each beam is 
an axial compression (or tension); no bending or twist- 
ing occurs. The "flow" of stress is parallel to the beams, 
as for the electric courant, and the geometry defined by 
the three conditions corresponds to the most homoge- 
neous distribution of constraints and currents through 
the whole structure. 

We must point out that condition c) is sufficient for 
having no bending in a structure for which conditions 
a) and b) are fulfilled, but not necessary: there do exist 
structures which do not satisfy condition c) and which de- 
form under compression without bending (e.g. see struc- 
tures of Fig. Nevertheless, the bulk modulus of such 
structures will be below the upper-bound value; condi- 
tion c) must be satisfied in order to have maximal bulk 
modulus. 

Furthermore, it is worth noticing that the three con- 
ditions are independent of the connectivity of the junc- 
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FIG. 3: Examples of cellular networks for wich condition c) is 
not satisfied and still deform by affine compression (no bend- 
ing or twisting of the beams). However their respective bulk 
modulus is strictly lower than the upper bound value SA. 



FIG. 4: In a real 3D foam, the edges (also called Plateau 
borders) meet in fourfold junction with equal angles, corre- 
sponding to the tetrahedral angle: arccos (—3) — 109, 5°. 



tions. As a consequence, there is an infinity of structures 
with maximal bulk modulus. Indeed, various examples 
of structures with maximal bulk modulus can be found in 
literature: we can cite the square, hexagonal, triangular, 
kagome networks H0 as 2D structures and the cubic 
and Kelvin networks as 3D structures. All 

these structures satisfy the necessary and sufficient con- 
ditions a), b), c), in agreement with the workpresented 
here. Numerical simulations of random 2D and 
3D 0] isotropic cellular materials have been also gener- 
ated. As expected, the bulk modulus of those materials 
is found to be always lower than the respective upper 
bound values. 

As concluding remarks, let us discuss the consequences 
of the criteria presented here on the mechanical proper- 
ties of real foams. Foam in the low-density limit is a 
particular cellular material: usually its preparation in- 
volves a continuous liquid phase that eventually solidifies. 
Therefore its structure is controlled by minimization of 
surface energy, leading to geometrical rules known as the 
Plateau's laws Q, which can be summarized as it follows: 

- edges in a 2D foam meet in threefold junctions with 
equal angles of 120°. 

- lamellae in a 3D foam meet in threefold lamella junc- 
tions (usually called Plateau borders) with equal angles 
of 120°, and Plateau borders meet in fourfold junctions 
with the tetrahedral angle: arccos (—5) — 109,5° (see 
Fig. EI). 

As a consequence, condition c) is always satisfied in 
a real foam. Usually, condition b) is nearly satisfied as 
well. The validity of condition a) is more delicate: while 
it is still possible to build 2D foams satisfying simulta- 
neously conditions a), b), c) and the Plateau's laws (e.g. 
the hexagonal honeycomb), this is no longer true for the 
3D case: no cell in a 3D foam is a simple polyhedron 
with straight edges, because a planar polygon cannot 
have all angles equal to the tetrahedral angle. Conse- 
quently, edges in a real foam must be curved, violating 
condition a) , and the bulk modulus and the average con- 
ductivity of a solid open-cell foam (i.e. a foam where the 



lamellae broke up during solidification) is always strictly 
lower that the corresponding bounds. While the conduc- 
tivity drop is not really significant when the edges are 
slightly curved [TlT ] , the bulk modulus value can be dra- 
matically decreased, because beams can easily bend or 
twist. We conclude that 3D cellular materials manufac- 
tured by some foaming process are probably not the most 
relevant for the design of high-bulk modulus/low-density 
structures. 

J. Robiche, J.-F. Sadoc and D. Weaire are deeply ac- 
knowledged for useful discussions and inspection of this 
work. 
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